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L'invention concerne les filtres numeriques realisant sur une 
fonclion frequentielle un traitement correspondanl a I'annulation d'une partie 
du signal temporel obtenu par transformde inverse d'une telle fonction 
frequentielle. 

Une telle fonction frequentielle peut elre la reponse frequentielle. 
c'est a dire la fonction de transfert d'un filtre. par exemple d'un filtre adaptatif. 
ou encore la transform6e de Fourier d'un signal temporel. 

De maniere connue. un filtre numerique realise des operations sur 
des donnees qui peuvent §tre soit des donnees directement disponibles a 
I'entree ou a la sortie du filtre. soit des donnees obtenues par une 
transformation telle que la transformation de Fourier, la transformation de 
Hartley, ou d'autres transformations. 

L'utilisation de ces transformations, frequemment realisees a I'aide 
d'algorithmes rapides, permet de reduire considerablement le nombre 
d'operations arithmeliques a realiser pour trailer les donnees. 

Les disposltifs et traitement numeriques sont souvent prevus pour 
realiser des operations de convolution ou de correlation. De nombreuses 
applications a base de traitement adaptatif du signal necessitent ainsi des 
operations de convolution (filtrage) et de correlation, qui sont realisees sur 
des vecteurs dont les donnees ou les coefficients sont variables, ces 
vecteurs etant souvent de grande taille (plusieurs centaines a plusieurs 

milliers de coefficients). 

Dans le domaine des telecommunications. I'annulation d'echo 
acoustique et electrique, le debnjitage de la parole, r§galisation de canal de 
transmission, et dans d'autres domaines, des applications comme la 
commande adaptative. font appel a de telles techniques de traitement 
adaptatif du signal. 

Le calcul exact de la convolution ou de la correlation necessite 
fr6quemment I'application d'une contrainte sur la taille du vecteur tempore! 
representant le filtre ou les donnees. c'est-a-dire une limitation du nombre 
d'echantillons non nuls de la transformee inverse du vecteur frequentiel a 
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utiliser. que ce vecteur frequentiel soit la reponse frequentielle d'un filtre ou la 
transformee d'un signal temporel a trailer. 

Cette conlrainte est classiquement realisee par application sur les 
donnees du domaine transforme. c'est-a-dire du domaine fr§quentiel. d'une 
5 transformation inverse, puis application d'un fenetrage annulant certaines 
composantes du vecteur temporel, et enfin application de la transformation 
directe pour obtenir le r6sultat exact dans le domaine transforme. 

Lorsqu'on souhaite realiser un filtre optimal § appliquer ^ un signal 
affecte d'une perturbation, pouvant fetre par exemple un bruit ou un echo, il 
1 0 est usuel de calculer la fonction de transfer! du filtre optimal a partir du signal 
affecte. sur un nombre de points frequentiels de la transformee du signal qui 
est inferieur au nombre de points total du signal. 

De plus, le vecteur frequentiel du filtre etant court, le calcul exact 
du signal de sortie du filtre, qui doit §tre realise par convolution lineaire. 
suppose souvent que la reponse impulsionnelle du filtre optimal soit 

prolong^e par des zeros. 

Cette operation de calcul d'un filtre a reponse frequentielle longue 
a partir d'un filtre a reponse frequentielle courte. en prolongeant la reponse 
impulsionnelle du filtre par des zeros, est habituellement effectuee. la encore, 
par passage intermediaire dans le domaine temps, prolongement du vecteur 
reponse impulsionnelle du filtre par des zeros, et transformation directe pour 
obtenir le resultat exact dans le domaine transforme. 

On comprend que ce passage dans le domaine temps pour 
appliquer une contrainte d'annulation, ou encore pour prolonger la reponse 
impulsionnelle d'un filtre optimal par des zeros, est couteux en calculs 
puisqu'il necessite deux transformations (une directe et une inverse). 

On va maintenant exposer les bases th§oriques montrant des 
raisons typiques pour lesquelles de telles contraintes sont souvent 
n^cessaires : 

L'utilite de la transformation de Fourier (TF) pour effectuer un 
filtrage ou pour calculer une correlation decoule du fait que les operations de 
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convolution et de correlation se traduisent par des simples multiplications 
dans le domaine transforme. Nous decrivons ici uniquement le cas de la 
convolution, le cas de la correlation est compl^tement analogue. 

Soit h(n) un filtre de longueur U et x(n) un signal d rentr6e de ce 
filtre. Alors on obtient a la sortie du filtre 



ou TF[..] et TFI[..] designent la transformation de Fourier directe et inverse et 

" designe I'operation de convolution. Dans la pratique, deux difficultes 

empechent une utilisation directe de cette equation : 
10 - pour calculer la transformation de Fourier du signal, on doit 

connaitre le signal entierement, c'est-a-dire tous les echantillons du passe et 

du futur, ce qui n'est pas realisable ; 

- dans de nombreuses applications comme I'annulation d'echo ou 

le filtrage adaptatif. le filtre h(n) evolue dans le temps. 
3^5 Pour resoudre ces problemes. I'homme de I'art decoupe le signal 

en blocs d'une longueur N et associe d chaque bloc b un filtre hb(n) de 

longueur U. Dans le cas general, ces blocs se chevauchent. Dans I'objectif 

de simplifier les notations, on suppose que les blocs sont juxtaposes, bien 

que ce ne soit pas une condition n6cessaire. 
.20 Soit 



Pour calculer la transformation de Fourier directe et inverse, on 
25 utilise une m§thode rapide connue sous le nom « Transformation de Fourier 
rapide » (TFRt pour une transformation directe sur T points. TFRIt pour la 
transformation inverse sur T points). 



>'(/!) = .v-(//) '//(//) = X''(')-v(//-') = 7/-7[/7-[.v(//)]V/ •[//(//)]] 




x{n + bN)siO^ n^N-^ 
0 ailleurs 



Le signal de sortie y(n) est donne par 



y(")= T.y>,i"-bN) avec 



yb(n)=Xb(n)*hb(n) 
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Le signal yb(n) est nul pour n<0 et n>T=N+U-1. ce qui implique 
que le nombre minimal de points pour calcuier les transformations 6gale T. 
En effet, on peut prouver que pour 0<n<T-1 on obtient : 
>..(/,) = 7™,[7Hi.k(/<)]- 7V--AV[/',.(")D = ■/•/•«/rl-\Vr(A-)- //,,,(A-)1 

Dans certains cas, on n'a pas besoin de toutes les valeurs de yb. 
mais uniquement de Ly points consecutifs. Dans ce cas. la taille des TFR doit 
satisfaire 

T>0.5*(N+U+Ly-1) 

Cette contrainte est appelee contrainte de convolution lineaire. 
Dans la pratique, on rencontre souvent le probl^me suivant : 
On s'est fixe N (la taille d"un bloc de donnees). T (la taille de la 
TFR) et Ly (longueur utile du resultat). On a calcule X(k) par TFR et le filtre 
H(k) par une methode arbitraire directement dans le domaine transforme sur 
T points. On dispose done de deux facteurs de la multiplication, mais la 
1 5 condition necessaire pour que le filtrage soit correct : 
U<2T-N-Ly+1 

n est pas remplie. La question est comment obtenir un filtre H{k) dans le 
domaine transforme qui soit une approximation de Hb.T(k). mais qui remplisse 
cette contrainte ? Dans la pratique connue. on calcule la transformee inverse 
de H. on annule une partie de ce signal et on recalcule un transformee 
directe de cette fonction prolongee pour remplacer H. 
On rencontre 6galement tr§s souvent un autre probleme, qui est le suivant : 

On s'est fix6 N (la taille d'un bloc de donn6es). T (la taille de la 
TFR) et Ly (longueur utile du resultat). On a calcule X(k) par TFR et le filtre 
H(k) par une methode arbitraire directement dans le domaine transforme sur 
L=T/D<2T-N-Ly+1 points. La question est comment obtenir un filtre H{k) 
dans le domaine transforme sur T points au lieu de L, qui soit une 
interpolation de H(k) et qui remplisse la contrainte Lh<2T-N-Ly+1 ? 
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Dans la pratique connue, on calcule la transformee inverse de H, 
on la prolonge par des zeros, et on recalcule une nouvelle fonction H a partir 
de celte transformee inverse tronquee. par une transformee directs. 

On a propose de supprimer ce passage par deux transformations 
successives. couteux en calculs. dans le cas d'une contrainte d'annulation. 
au moyen d'un traitement (filtrage) directement dans le domaine transforme 
sur les donnees du domaine frequentiel. 

Dans le domaine de I'annulation de Techo acoustique, et plus 
precisement dans le domaine des filtres adaptatif par blocs avec utilisation 
de la transformee de Fourier, appeles frequemment « filtres adaptatifs dans 
le domaine frequences » {« Frequency Domain Adaptive Filter », note FDAF 
), Clark et ColL [1] decrivent un formalisme general ou apparait la necessite 
d'appliquer une contrainte sur la longueur de la reponse impulsionnelle du 
filtre adaptatif pour obtenir un algorithme exact, et proposent de mettre en 
application une convolution avec la reponse frequentielte du filtre. qui soil 
equivalents a une multiplication par une fenetre dans le domaine temporel. 
Le nombre de coefficients de la convolution proposee est tres eleve. 
maintenant une grande complexity de calcul. 

Sommen et Coll. [5] decrivent une autre variante de FDAF ou la 
contrainte temporeile est realises par convolution dans le domaine 
frequentiel : les auteurs proposent d'appliquer sur les donnees frequentielles 
representant la reponse du filtre un lissage particulier a 3 coefficients 
equivalent a une ponderation de la reponse temporeile par une fenetre en 
cosinus. 

Le maximum de cette fen§tre est suppose correspondre aux 
coefficients dominants de cette reponse. ce qui est une hypothese tres 
restrictive. La justification donnee par les auteurs a leur methode est 
d'ailleurs la reduction du bruit d'adaptation liee a la ponderation des 
coefficients, plutot que la recherche d'une approximation de la solution 
exacte d'application d'une contrainte. 



On citera egalement a titre de references. Mansour et Coll. [2], qui 
ont propose une version simplifiee du FDAF sans application de contrainte. 
ce qui se traduit par une reduction importante de complexity, mais aussi par 
des performances inf^rieures a celles des algorithmes avec contraintes. 
5 Prado et Coll. [3] ont egalement decrit une version generalisee du 

FDAF appelee GMDFcx (« Generalized Multi-Delay Filter » avec facteur de 
surechantillonnage a), qui permet de Iraiter des reponses impulsionnelles 
longues avec des blocs de donnees plus courts en segmentant la reponse 
impulsionnelle en petits blocs, ce qui reduit le retard de traitement dans 

10 Talgorithme. Get algorithme utilise la contrainte citee precedemment. 

Mc Laughlin et Coll. [4] decrivent une forme modifiee de 
I'algorithme FDAF reprenant certains aspects du GMDF, notamment la 
segmentation de la reponse impulsionnelle en petits blocs, et proposent une 
technique d'application de la contrainte sur les differents blocs evoluant au 

15 cours du temps suivant un schema bien defini (« scheduler »), qui permet de 
limiter la complexite tout en preservant partiellement le benefice de la 
contrainte. 

Le but principal de ['invention est de proposer un filtre numerique 
effectuant, sur un vecteur frequentiel, une convolution qui equivaut a un 
2 0 fenetrage du signal temporel correspondant a ce vecteur frequentiel, qui 
produise une precision de calcul satisfaisante tout en ne necessitant qu'une 
quantite de calculs particulierement faible. 

Ce but est atteint selon I'invention grace a un dispositif de 
traitement de donnees numeriques comportant un module realisant sur un 

2 5 vecteur de donnees du domaine frequentiel Z(k), ou k varie de 0 a N-1, une 

convolution avec une fonction U, convolution qui correspond a une 
annulation dans le domaine temporel d'echantillons de la transformee 
inverse de Z(k), caracterise en ce que la fonction U est de la forme : 

U{k) = sine ~ B ^ ^ P(k) , ou koest un entier constant et P(k) est 

V 2 ^ 

3 0 une fenetre de ponderation symetrique autour de ko. 



Un tel fillre presente un nombre de coefficients qui peut etre ajuste 
suivant rapproximation desiree du resultat, ce nombre pouvant etre tres 
faible. Un tel filtre ne necessite. pour etre mis en oeuvre. qu'un nombre 
d'operations arithmetiques considerablement reduit par rapport aux solutions 
classiques. 

D'autres caracteristiques. buts et avantages de I'invention 
apparaitront a la lecture de la description qui va suivre. faite en reference aux 
figures annexees sur lesquelles : 

- les figures 1 a 3 representent respectivement la reponse 
frequentielle en module, la reponse frequentielle en phase, et 
la reponse impulsionnelle d'un filtre, 

- les figures 4 a 6 representent respectivement la reponse 
frequentielle en module, la reponse frequentielle en phase, et 
la reponse impulsionnelle du filtre des figures 1 a 3 sur lequel 
une meme contrainte d'annulation a ete appliquee a I'aide d*un 
dispositif classique a deux transformees, 

- les figures 7 a 9 representent respectivement la reponse 
frequentielle en module, la reponse frequentielle en phase, et 
la reponse impulsionnelle du filtre des figures 1 a 3 sur lequel 
une contrainte d'annulation a ete appliquee avec un dispositif 
selon rinvention, 

- les figures 10 a 12 representent respectivement en module 
frequentiel. en phase frequentielle, et en reponse 
impulsionnelle I'erreur obtenue avec le dispositif selon 
I'invention par rapport au dispositif classique, 

- la figure 13 represente I'amplitude de la fonction de transfer! 
d'un filtre et quelques points interpoles avec un dispositif selon 
I'invention. 

- la figure 14 represente la reponse impulsionnelle du filtre de la 
figure 13 non interpole, prolongee par des zeros. 



la figure 15 represente la reponse impulsionnelle du filtre de la 
figure 13 interpole. 

la figure 16 represente I'erreur en reponse impulsionnelle entre 
le filtre de la figure 14 prolonge par des zeros dans le domaine 
temps et le filtre interpole de la figure 13, 

les figures 17 a 21 representent cinq variantes de dispositifs a 
filtre adaptatif conformes a rinvention, 

- la figure 22 represente schematiquement un annuleur d'echo 
de type connu, 

- la figure 23 represente de maniere fonctionnelle detaillee un 
annuleur d'echo connu, 

la figure 24 represente de maniere fonctionnelle detaillee un 
annuleur d'echo salon rinvention. 

- la figure 25 represente un dispositif a annulation d'echo et de 
bruit, 

- la figure 26 represente un dispositif ayant le meme role que 
celui de la figure 25, et realise selon rinvention. 

On decrira maintenant un premier dispositif de filtration d'un son 
conforme a rinvention en reference a la figure 17. 

Ce dispositif re^oit en entree un signal temporel discret constitue 
de plusieurs blocs consistant chacun en 1024 echantillons dans I'exemple 
non limitatif decht ici. Ce signal temporel est un signal sonore ou acoustique, 
tel qu'un son representatif de la parole d'un utilisateur parlant dans un micro. 

Sur la figure 20. un de ces blocs temporels est note S1(0,....N-1). 
ou N est egal a 1024. 

Ce signal SI est regu par un premier module Ml qui realise sur ce 
signal une transformee de Fourier discrete rapide d'ordre N, notee FFTn. 
pour fournir un signal frequentiel de N 6chantillons X(0,.,.,N-1). Ce signal 
frequentiel est done lui aussi un signal representatif d'un son. Ce signal est 
done ici un signal acoustique. 



A partir de ce vecteur frequentiel X de taille N. un filtre adaptatif H 

calcule un vecteur H(0 N-1) dont les N coefficients dependent du vecteur 

X. Chaque composante du vecteur H est calculee gr^ce a un algorithme 
faisant intervenir des composantes du vecteur X. Cette dependance est 
indiquee sur la figure 17 par une double fidche. 

De manidre connue. un tel filtre adaptatif H est destine a etre 
applique au signal X ^ partir duquel il s*adapte. par exemple pour annuler un 
echo detecte sur ce signal X. 

Toutefois, le filtre H a N coefficients ne peut pas etre applique a X 
directement tel qu'il est obtenu en sortie de H. Comme explique 
precedemment. le vecteur H doit etre traite de fafon a ce que sa reponse 
impulsionnelle presente un nombre de coefficients non nul qui soit inferieur a 
une valeur N' qui est inferieure a N. En d'autres termes. le vecteur H doit 
repondre a une contrainte de convolution lineaire avant d'etre convolue avec 
le vecteur X. 

Pour obtenir un vecteur H* repondant a une telle contrainte a partir 
du vecteur H, il est connu de placer en sortie du filtre adaptatif H un 
ensemble de trois modules. Tun effectuant une transformee inverse de H 
dans le domalne temps, un deuxieme module multipliant la reponse 
impulsionnelle du filtre H par une fenetre rectangulaire qui garde inchangee 
une partie de la reponse impulsionnelle et annule une partie de queue de 
cette reponse impulsionnelle. et un troisieme module effectuant une 
transformee directe pour obtenir la reponse frequentielle du filtre 
correspondent a la reponse impulsionnelle tronquee. 

Le dispositif de la figure 17 presente entre la sortie du filtre 
adaptatif H et un module de convolution avec le signal X, un module M3 
realisant sur le vecteur frequentiel H(0.....N-1) un traitemenl dans le domaine 
frequentiel. sans passage par le domaine temporel, pour fournir un vecteur 
frequentiel H' a N coefficients dont la reponse impulsionnelle n'a pour valeurs 
non nulles que ses N' premiers coefficients, ces N' premiers coefficients 
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etant les memes que les N' premiers coefficients de la reponse 
impulsionnelle de H. 

Pour effectuer un tel trailement, ce module M3 met en oeuvre une 
convolution de H(0,...,N-1) avec une fonction U(k) qui est donn6e par 

Uik) - -sincf-^^^^-i^ ' pour k allant de 0 a Lu-1. 

et U(k) etant nul en dehors des valeurs de k allant de 0 a Lu-1 . La fonction 

, ^ sin(;r.x) , - 

sine de la formule ci dessus est donnee par : sm c(x) = . La fonction 

TT.X 

sine prend la valeur 1 pour x=0. 

Dans rexemple particulier decrit ici. Lu est choisi egal a 7. de sorte que Lu-1 
est pair, et de sorte que (Lu-1)/4 est un multiple de 1/2. 

La fonction U(k) peut encore s'ecrire 

Lf{/c) = ^.sinc(^^^j.e Kaiser (k ^ fl) , 

Une telle fonction U(k) a pour transformee temporelle une fenetre 
sensiblement reclangulaire. prenant pour valeur environ 1 pour k allant de 0 
a 512, et pour valeur 0 pour k allant de 512 a 1024. 

Les inventeurs sont arrives a une telle forme de fonction en 

constatant que la fonction sin c - a pour transformee discrete inverse dans 

le domaine temporel, une fonction prenant une valeur sensiblement 
constante et proche de 2 pour k compris entre 0 et N/4 ainsi que pour k 
compris entre 3N/4 et N, tandis qu'elle prend des valeurs proches de zero 
pour k compris entre N/4 et 3N/4. 

Pour obtenir une fonction frequentielle U(k) correspondant a une 
fenetre temporelle prenant sensiblement la valeur 1 pour k=0,...,N/2, et la 

valeur 0 pour k=N/2 N-1, les inventeurs ont utilise la propriete selon 

laquelle une multiplication dans le domaine frequentiel par une exponentielle 



complexee ' ^ . equivaut a un decalage de no de la fonction temporelle 
correspondante, obtenue par transformee de Fourier inverse. 



produisant un decalage circulaire d'un quart de fenetre de la transformee 
discrete inverse de ta fonction sin c| - 1. permettant done d'obtenir la fenetre 



souhaitee. Cette fenetre correspond done a I'annulation d'une deuxieme 
moitie du signal temporel. et a la conservation telle quelle de la premiere 
moitie du signal temporel. Plus generalement on pourra utiliser un decalage 
different de celui d*un quart de fenetre propose ici. On utilise alors une 



choisie pour obtenir le decalage souhaite. Ainsi, une valeur de a egale a 2 
correspond a une fenetre temporelle qui est nulle sur le premier et le 
quatrieme quart du support temporel. et egale a 1 sur une partie centrale du 
support, de longueur egal a la moitie de celle du support. L'invention ne se 
limite done pas a une valeur de a egale a 1 comme on le decrit dans les 
presents exemples particuliers de realisation. Un autre exemple preferentiel 
de realisation de I'invention est celui ou a = -1. correspondant a une 
annulation de la premiere partie de la fenetre. L'homme de Tart adaptera la 
valeur de a selpn la partie qu'il souhaite annuler, en tenant compte que le 
decalage se fait de maniere circulaire. 

Dans cette formule. la fonction U est multipliee par une fenetre de 
ponderation classique. ici la fenetre de kaiser avec p=1 .5. 

La fenetre de ponderation a pour role connu en lui-meme de 
reduire un pic d'amplitude d'erreur sur la fenetre transformee d'un signal 
tronque, ici du signal 

U{k) = sin c[ ^ " le ^ ^ tronque sur I'intervalle 0 ^ /c < (La - 1). 



Dans le cas present, les invenleurs ont utilise la valeur hqIN =1/4. 





[kaiser {k ^ fi) , ou a est une constante 



De maniere preferentielle, la fenetre de ponderation est une 
fenetre symetrique par rapport a une valeur ko de k pour laquelle la fonction 

U{k) sine 



V 2 J 



presente son maximum en module (en valeur 



U{k) ^ sill cf j.e ^ et de cetle fenetre de ponderation presentent 



absolue). De cette fa^on. les coefficients du produit de 
2 j 

les memes proprietes que les coefficients de L/(/c) = sin d — ^ J.e 

qui seront decrites ci-apres. et qui procurent des avantages importants en 
termes de quantite de calcul et de taille memoire. 

L'homme de Tart peut done remplacer la fenetre de kaiser ici 
proposee par toute autre fenetre de ponderation connue. telle qu'une fenetre 
de Manning. 

Dans I'exemple de realisation present, avec Lu=7. on obtient pour 
les coefficients de U suivants : 

k 0 1 2 3 4 5 6 

U(k) j 0.0644 0 j 0.3027 .5 -j 0.3027 0 -j 0,0644 

On constate que le filtre presente une symetrie hermitienne par 

rapport a la valeur ko de k pour laquelle - ' + y = 0 , c'est a dire 

k^ = ^^^^ . c*est-^-dire ko = 3. 
0 2 

Grace a cette symetrie hermitienne en elle-meme. seule la moitie 
des coefficients du filtre U sont a calculer. De plus. grSce a cette symetrie 
hermitienne en elle-meme, ['operation de convolution avec U pourra etre 
effectuee avec une reduction de la charge de calcul qui est proche de 50%, 
comme on le montrera par la suite. 

On constate en outre que les coefficients du filtre sont 
alternativement reels purs et imaginaires purs, Cette propriete vient du fait 
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que I'argument du facteur e * - est une fonction lineaire variant par pas 

de -i - . et a pour valeur ^ k=0 un multiple de j ^ . 
2 - 

Cetle propriete vient du fait que le filtre est un filtre demi-bande. 
Toutefois. une symelrie hermitienne peut egalement etre obtenue avec un 
autre filtre qu'un filtre demi-bande, par exemple un filtre correspondant a une 
fenetre annulant un quart du signal temporel. 

On constate encore qu'un coefficient sur deux du filtre U est nul. 
sauf le coefficient pour k=3. En effet. la fonction sinc(x) s'annule pour tout x 
entier sauf pour 0. Cetle propriete decoule egalement du fait que U est un 
filtre demi-bande. Ainsi. dans le cas des coefficients indiques precedemment, 
un coefficient sur sept est nul. 

Ces trois proprietes permettent done une reduction importante de 

la charge de calcul. 

La convolution de U(k) avec H(k) realisee par le module M3 pour 

fournir le vecteur H' s'exprime sous la forme : 

H'(k+3) = H(k).U(6)+ H(k+1).U(5) + H(k+2).U(4) + H(k+3).U(3)+ 
H(k+4).U(2) + H(k+5).U(1)+H(k+6).U(0) 

On calcule done un coefficient d'indice k' de H' en convoluant U et 
H de telle sorte que le coefficient de H de meme indice k* est multiplie par la 
valeur de U qui est maximaie en module. Les autres multiplications de cette 
convolution sont realisees en augmentant I'indice de U lorsque I'indice de H 
est diminue, comme pour toute convolution. 

En d'autres termes. la fonction H' est egale a la convolution de H 

avec une fonction de forme U{k) = ^ . sin cf ^^^-2- J.e ' ^ ou ko est nul. 

Plus generalement. une fonction U de la forme 
t;(/c) = ^sinc|^S + ^ + |je'"f^"''^lp(/c). ou la constante B est un entier 
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choisi de maniere arbitraire et ou P est une fenetre de ponderation. peut etre 
utilisee. 

En effet. la constante B produit un decalage de la fonction U(k) qui 
sera facilement pris en compte par Thomme du metier dans le calcul de la 
convolution, en tenant compte du fait que pour obtenir le coefficient d'indice 
k' donne de H\ la convolution doit Stre telle que le coefficient de H de meme 
indice k' doit 6tre multiplie par le coefficient de U de module maximal. 

Le fait qu'une fonction F2 soil egale a une convolution entre F1 et 
Uo. signifie que : 

F2(k)= F1(k).Uo(0) + F1(k-1).Uo(1) +F1(k+1).Uo(-1) + F1 (k-2).Uo(2) + 
F1(k+2).Uo(-2) + ... 

Dans le cas ou U (k) = Uo(k+B), on peut ecrire : 

F2(k)= F1(k).U(-B) + F1(k-1).U(-B+1) +F1(k+1).U(-B-1) + ... 
Etdonc F2(k+B)= F1(k+B).U(-B) + F1{k+B-1 ).U(-B+1 ) +F1(k+B+1).U(-B-1) 

+ ... 

II est done sans consequence sur le resultat d'une convolution de 
decaler Tune des fonctions convoluees d'un nombre entier d'echantillons. si 
le signal obtenu par convolution est lui-meme decale de maniere appropriee. 

Le module M3 peut done se decomposer en un premier ensemble 
realisant les multiplications et les additions, et un second ensemble realisant 
le decalage. 

Dans le cas present et conformement a invention, on effectue 
une convolution entre le signal frequentiel d'entree H et la fonction U de 
fa^on a ce que, dans le calcul de convolution donnant H' a I'indice k. H(k) est 
multipliee par la valeur atteinte par U lorsque le sine de U est maximal, c'est 
a dire que H(k) est multipliee par U(k+A) ou A est tel que ['argument du sinus 
cardinal est nul. 

Dans cette convolution, H(k+1) est done alors multipliee par U(k+A-1), H(k-1) 
par U(k+A+1), H(k+2) par U(k+A-2). H(k-2) par U(k+A+2), et plus 
generalement H(k+p) est multiplie par U(k+A-p). 
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Les inventeurs ont constate qu'en realisanl une telle convolution 
avec une telle fonction U. on realise Tapplication d'une contrainte 
d'annulation de la moitie du signal temporel correspondant a H. avec une 
precision tres satisfaisante. et pour une quantite de calcul tres faible. sans 
avoir a realiser une transformee inverse dans le domaine temporel . 
Dans le cas du filtre U de Texemple de realisation de la figure 17, on a. 
comme expose prec6demment. 

U(l)=U{5)=0. U(0) =-U(6). U(1) =-U(5). et U(2) =-U(4). 
Done H'(k+3)=U(0).( H(k+6)-H(k))+U(2)( H(k+4)-H(k+2))+ U(3)(H(k+3)). c'est 
a dire H'(k+3)=U(0).A(k)+U{2)B(k)+ U(3)(H(k+3)) avec A(k) = H(k+6)-H(k) et 
B(k) = H(k+4)-H(k+2). 

Ainsi, au lieu de 7 multiplications complexes et 6 additions 
complexes, c'est-a-dire 28 multiplications reelles et 26 additions reelles. on 
realise la convolution avec 6 multiplications reelles et 8 additions reelles. 

Dans le cas d'un dispositif different de celui de la figure 17 mais 
presentant un module M3 similaire. si le signal H est obtenu par transformee 
d'un signal temporel reel, alors H presente une symetrie hermitienne pouvant 
elle aussi etre mise a profit et permettant encore de diviser la charge de 
calcul par 2. 

Les valeurs H(k) pour k<0 et k>N-1 peuvent etre obtenues par un 
adressage modulo N (-1 correspond a N-1, -2 a N-2, etc.). 

On a represents sur les figures 1 et 2 le signal H(k) en amplitude 
et en phase sur I'ensemble de ses 1024 points, et sur la figure 3, la reponse 
impulsionnelle h(i) sur 1024 points, correspondant a H(k). 

Les figures 7 et 8 representent le signal frequentiel H' obtenu 
apres convolution de H et U. respectivement en amplitude et en phase. La 
figure 9 represente la reponse impulsionnelle h' correspondant a ce signal. 

On constate que la reponse impulsionnelle h' est annul6e sur les 
512 derniers points et elle reste inchangee sur ses 512 premiers points. 

On compare maintenant ces resultats obtenus avec le dispositif de 
rinvention, avec les resultats obtenus avec un dispositif classique mettant en 
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oeuvre une transformee de Fourier inverse, une annulation de la deuxieme 
portion de la reponse impulsionnelle, et une transformee de Fourier de ce 
signal temporel tronque. 

Les figures 4 et 5 repr6sentent ainsi le signal frdquentiel obtenu 
5 apr^s ces trois etapes de traitement classiques. respeclivement en amplitude 
et en phase. La figure 12 represente la reponse impulsionnelle de ce signal. 

On a represente sur les figures 10 et 11 I'erreur respeclivement en 
amplitude et en phase entre les signaux frequentiels obtenus avec ce 
dispositif classique et les signaux frequentiels obtenus avec le dispositif 
10 selon invention. 

La figure 12 represente I'erreur sur la reponse impulsionnelle. 

L'homme de I'art observera que I'erreur obtenue est tres faible. 
compatible avec les performances requises notamment dans les applications 
de traitement de parole telles que I'annuiation d'echo acoustique. 
15 On decrira maintenant un deuxieme dispositif selon I'invention, en 

reference a la figure 18. Ce dispositif re?oit en entree un signal temporel SI 
a 2N echantillons. 2N etant ici egal a 1024. Un premier module Ml realise 
une transformee de Fourier discrete a 2N points de SI pour fournir un 
vecteur frequentiel X a 2N points, d'indices allant de 0 a 2N-1 , 
2 0 Un module M5 place en sortie de Ml realise une decimation de 

X(0,....2N-1). de sorte qu'il ne garde qu'un coefficient sur deux de X. pour, 
constituer un vecteur frequentiel X'(0,...,N-1). 

En sortie de M5 est place un filtre adaptatif dependant de X' et 
fournissant un vecteur H(0„..N-1) a N coefficients. Le vecteur H ne 

2 5 comportant que N coefficients, il ne peut pas etre convolue avec X(0,....2N- 

1). II est done necessaire de construire a partir de H un vecteur H'(0....,2N-1) 
a 2N coefficients dont la transformee inverse ait ses N premiers 6chantillons 
identiques a ceux de la reponse impulsionnelle de H, et dont les N derniers 
echantillons soient nuls. 

3 0 On place pour cela un module M3 entre le filtre H et un module de 

convolution avec la fonction X(0,...,2N-1). 
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Le module M3 re^oit en entree le signal frequentiel H a N=512 
points, correspondant a une reponse impulsionnelle de 512 points. 

Partant de ces N=512 points frequentiels, on cherche a obtenir 
une fonction H' ^ 2N points. 

Le module M3 se decompose en deux modules M3* et M3". Le 
premier module M3' insere un zero apres chacun des N points du signal H 
pour obtenir un signal Ho(0.... .2N-1 ) a 2N = 1024 points dont un coefficient sur 
deux est nul, et plus precisement dont les coefficients d'indices impairs sont 
nuls. 

Le vecteur frequentiel Ho ainsi obtenu correspond a une reponse 
impulsionnelle a 2N=1024 points dont les N premiers points sont les memes 
que ceux de la reponse impulsionnelle de H a un facteur pres. 

La reponse impulsionnelle de H reste done inchangee par 
rintroduction de zeros entre les echantillons initiaux, si les echantillons 
initiaux se retrouvent dans le signal complete a un echantillon sur deux en 
partant du premier indice du signal complete, c'est a dire a partir de I'indice 0 
si le signal complete est indice a partir de I'indice 0. Les echantillons de H se 
retrouvent done aux indices pairs de H'. Une telle insertion ne modifie pas la 
premiere partie a 512 points de la reponse impulsionnelle. 

Un second module M3" realise ensuite un traitement du vecteur 
frequentiel complete Ho. qui est similaire au traitement decrit prec6demment 
en reference a la figure 17, et qui a pour effet d'annuler les N derniers points 
de la reponse impulsionnelle a 2N points de Ho sans modifier les N premiers 
points, et cela a partir du domaine transforme sans revenir dans le domaine 
temporel. 

On utilise ici un filtre de longueur 15 dont les coefficients se 
calculent d'une maniere analogue a ceux du premier exemple de realisation. 

Dans le cas present, une valeur sur deux du signal frequentiel Hoa 
convoluer est nulle. Ainsi. seulement un coefficient sur 2 du filtre doit etre 
multiplie avec une valeur non nulle de Ho(k). 



De plus, le fait d'annuler la seconde partie de la reponse 
impuisionnelle et de conserver inchangee sa premiere moitie ne modifie pas 
les 512 coefficients initiaux de la fonction frequentielle entre lesquels ont et6 
introduits ^ chaque fois un zero de la fa?on prec6demment d6taill6e. 

Ainsi. la convolution de Ho avec une fonction U telle que celle 
proposes pr6c6demment. correspondant ^ une fendtre temporelle demi- 
bande. ne modifie pas les 512 valeurs non nulles de Ho. c'est ^ dire les 
valeurs Ho(k) pour les valeurs paires de k allant de 0 ^ 1022. 

Au bilan. seules les valeurs nulles de Ho(k) sont a recalculer par 

convolution. 

II n'est done necessaire que de calculer un coefficient sur deux du 
filtre U. Avec un filtre U de longueur 15 et atteignant son maximum en 
module, c'est a dire son point de symetrie hermitienne, pour ko=7. c'est a dire 

avec ^7(A) = sin/-^^^ + -V"'^^''(*'''^''") ®^ Uik) = rma "^ - ''K pour 
\^ A 2J 

Lu = 1 5. on obtient les coefficients suivants : 

K 0 2 4 6 7 8 10 12 14 

U(k) 0.0552j 0.1001 j 0.1952j 0.6308j 1 -0.6308j -0.1952j -O.IOOIj -0.0552) 

On choisit de preference Lu impair pour obtenir autant de 
coefficients de part et d'autre du point de symetrie de U. 

Un coefficient sur deux etant inutile, les coefficients a calculer sont 
tous purement imaginaires et presentent bien une symetrie hermitienne par 
rapport a ko=7. 

Comme dans le cas precedent, les calculs sont avantageusement effectues 

de la fa^on suivante : 

En posant A(k) = Z(k+1 4)-Z(k) 

B(k) = Z(k+12)-Z(k+2) 

C(k) = Z(k+10)-Z(k+4) 

D(k) = Z(k+8) - Z(k+6) 

et avec /<„ = ^^^^ = 7 on obtient le resultat interpole : 
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^H{k)si kest pair 

^^^^ " \U{0)A{k - 7) + U{2)B{k - 7) + U{4)C{k - 7) + U{6)D{k - 7)si k impair 

La complexite de calcul est done de 8 multiplications et de 14 
additions reelles pour chacun des points calcules. 

Bien entendu, dans un dispositif different de celui de la figure 18 
mais comportant un module M3 similaire. le vecteur H ^ interpoler peut §tre 
la transformee du signal temporel et ce signal temporel peut etre reel. Dans 
un tel cas. la complexite est encore reduite d'un facteur de 2 par la symetrie 
hermitienne de H. 

Nous montrons maintenanl les resultats correspondant a cet 

exemple. 

Sur les figures annexees sonfrepresentees : 

- Sur la figure 13. on a represente le module de la fonction de 
transfert du filtre H sur quelques points ainsi que quelques 
points interpoles pour former le vecteur H' ; 

- sur la figure 14. la reponse impulsionnelle d'un filtre H" 
obtenue a partir de celle de H par la methode connue 
consistent en deux transformees de Fourier rapides entre 
lesquelles est effectue un prolongement de la reponse 
impulsionnelle de H par 512 zeros, 

- sur la figure 15, la reponse impulsionnelle du filtre H' obtenu en 
sortie du module M3 du dispositif de la figure 18, 

- sur la figure 16, 1'erreur entre la reponse impulsionnelle du filtre 
H" et la reponse impulsionnelle du filtre H'. 

L'homme de Tart observera que I'erreur obtenue est, la aussi. 
compatible avec les performances requises notamment dans les applications 
de transmission de parole telles que le d6bruitage par filtrage optimal. 

Dans le dispositif de {'invention, on remplace done la suite 
d'operations complexe : transformation inverse, fenetrage ou prolongement 
par des zeros, puis transformation directe. par un filtrage approprie. Ce 
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filtrage de faible complexite est applique directement a des donnees issues 
de la transformation, qui est frequemment la transformation de Fourier. 

Dans une premidre variante, le dispositif permet d'approximer. 
avec un degr6 de precision adapte ^ I'application, un fen§trage temporel 
mettant a zero M echantillons consecutifs d'un bloc de T echantillons (M<T). 
Le filtrage appliqu6 aux donn6es dans le domaine transforme peut etre 
considere comma I'^quivalent d'un filtrage passe-bas, c'est-a-dire d'un 
lissage. 

Dans une seconde variante, le dispositif permet d'approximer. 
avec un degre de precision adapte a I'application, le prolongement d'un bloc 
de N echantillons par M zeros. Le filtrage applique aux donnees dans le 
domaine transforme peut alors etre considere comme une interpolation 
permettant de construire les M points manquants pour prolonger les donnees 
frequentielles a N+M points. 

Les figures 19, 20 et 21 montrent d'autres exemples de dispositifs 

conformes a invention. 

Chacun des dispositifs des figures 19. 20 et 21 presente un filtre H 
adaptatif a coefficients variables, a N ou 2N coefficients. Chacun de ces 
dispositifs presente un module M1 recevant en entree un signal temporel SI 
et realisant une transformee de Fourier de ce signal. 

Le premier module Ml du dispositif de la figure 19 regoit en entree 
un bloc temporel SI a N Echantillons. Ce module M1 realise une transformee 
de Fourier discrete a N points sur ce signal SI et fournit ainsi la transformee 
de Fourier X de ce bloc de N donnees. 

Dans le cas du dispositif de la figure 20, le module Ml complete le 
signal SI par des zeros pour obtenir un signal a 2N coefficients, puis realise 
la transformee de Fourier discrete a 2N coefficients de ce signal complete. 

Dans le cas de ces deux dispositifs, la fonction de transfer! du 
filtre H est calculee a partir de la transformee de Fourier X en sortie du 
module M1 et s'adapte en fonction de cette transformee de Fourier X. Cette 
dependence est symbolisee par une double fleche verticale. 
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Dans le cas de la figure 19. le dispositif presente en sortie du 
module M1. un module M2 qui realise les deux operations d6crites 
precedemment d'introduction d'un zero entre chaque couple de valeurs du 
signal X ^ N coefficients, pour obtenir un signal XX0.....2N-1). puis realise sur 
le signal X* une convolution avec un filtre U tel que ceux proposes 
precedemment. 

Ce module M2 realise done une interpolation pour etendre 
directement la taille de la transformee de Fourier a 2N points frequentiels a 
partir du signal X a N points. X* ayant pour reponse impulsionnelle la reponse 
impulsionnelle de X prolongee par des zeros. 

Le dispositif de la figure 19 presente egalement un module M3 
realisant sur la fonction de transfert H du filtre a N points la meme operation 
d'interpolation, c'est a dire Tintroduction d'un zero entre chaque couple de 
valeurs de H, puis une convolution avec une fonction U telle que celles 
proposees precedemment. correspondant a une annulation d'une deuxieme 
partie de la reponse impulsionnelle du filtre a 2N points obtenu par 
introduction des zeros. 

Ainsi. sur la figure 19. on obtient en sortie du module M2 un signal 
frequentiel representatif d'un signal temporel a 2N points dont les N premiers 
sont ceux du signal d'entree S1. On obtient en sortie du module M3 un signal 
frequentiel S3 a 2N points dont la reponse impulsionnelle a 2N points 
presente pour ses N premiers points les memes valeurs que la reponse 
impulsionnelle du filtre H, et pour ses N derniers points des valeurs nulles. 

Ces deux signaux frequentiels a 2N points S2 et S3 sont ensuite 
multiplies au niveau d'un module multiplicateur M4, qui peut §galement etre 
un module realisant une correlation ou une convolution entre S2 et S3. 

Le dispositif de la figure 20 pr6sente une structure similaire a celui 
de la figure 18. ^ la difference que le filtre adaptatif H presente 2N 
coefficients, et que le module M3 annule les N derniers coefficients de la 
reponse impulsionnelle du filtre a partir du domaine frequentiel. Le lissage 
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permet ainsi de realiser I'equivalent de la contrainte de convolution lineaire 
en limitant la longueur de la reponse impulsionnelle du filtre a N points. 

La reponse frequentielle H' ainsi obtenue est multipli6e au vecteur 
X au niveau d'un module muitiplicateur de sortie M4. 
5 La figure 21 represents un dispositif de m§me structure que celui 

de la figure 18, a la seule difference que le module Ml fournit une 
transformee frequentielle ^ N 6chantillons, et que le d6cimeur effectue une 
decimation de la transformee de Fourier X a N echantillons du bloc de N 
donnees en ne conservant qu'un point sur D (D etant un nombre entier 
1 0 superieur ou egal a 2). 

On obtient ainsi un vecteur X* a N/D points, a partir duquel on 
effectue I'adaptation de la fonction de transfert du filtre H, cette fonction de 
transfert H ayant L=N/D points. De la meme maniere que pour le dispositif de 
la figure 18, on reconstitue par interpolation une fonction de transfert H' sur N 
15 points au niveau d'un module M3 conforme a ceux exposes en detail 
precedemment, avant d'effectuer une convolution lineaire entre H' et X au 
niveau d'un module muitiplicateur M5. 

Au sein de M3. un premier module M3' introduit apres chaque 
echantillon de H un nombre D de zeros. Le signal frequentiel Ho ainsi obtenu 
20 correspond a une reponse impulsionnelle a D.L points dont les L premiers 
points sont les memes que ceux de la reponse impulsionnelle de H. La 
reponse impulsionnelle de H reste en effet inchangee par Tintroduction de 
zeros entre les echantillons initiaux. 

Un deuxieme sous-module M3" realise ensuite un traitement du 

2 5 signal frequentiel Ho a D.L points obtenu, qui a pour effet d'annuler les (D- 

1),L derniers points de la reponse impulsionnelle a D.L points de Ho . sans 
modifier les L premiers points de cette reponse impulsionnelle. On obtient 
ainsi le vecteur frequentiel H' souhaite a N=D.L points. 

Un telle decimation a I'avantage de reduire sensiblement la 

3 0 complexite de calcul pour adapter le filtre H ; elle est applicable dans les cas 
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OU la fonclion de transfert du filtre est assez reguliere pour etre decimee sans 
perte critique d'information. 

On notera que. grSce a invention, des filtres a tres petit nombre 
de coefficients peuvent §tre utilises. 

On decrira maintenant un dispositif d'annulation de I'echo 
acoustique et un dispositif de reduction d'une perturbation (bruit et echo) 
mettant en oeuvre un dispositif de I'invention. 

II est bien connu que dans certaines applications de transmission 
de la parole, telles que la telephonie mains-libres et la teleconference, des 
dispositifs appropries doivent etre utilises pour combattre I'echo acoustique. 
qui serait sinon une source de gene considerable dans les conversations. 

Dans I'etat de I'art. ces dispositifs mettent en oeuvre comme 
represents a la figure 22. un filtre adaptatif FA dont le role est d'identifier et 
de modeliser la reponse impulsionnelle du chemin de couplage acoustique 
entre un haut-parleur HP et un microphone MIC du terminal telephonique 
mains llbres ou du syst6me de teleconference. 

Cette partie du dispositif representee figure 22. est appeiee 
annuleur d'echo acoustique. 

Les realisations pratiques de ces dispositifs peuvent prendre des 
formes tres diverses. par exemple celles expos6es dans le document [8]. 

Les coefficients du filtre adaptatif FA sonl ajustes au cours du 
temps par un algorithme approprie. qui exploite un signal de parole reiju et 
une estimation d'erreur presente sur le signal de parole emise a partir du 
micro MIC et apres application sur ce signal du filtre FA. Un algorithme 
dassique pour I'adaptation des coefficients du filtre FA est connu sous le 
nom de gradient stochastique ou LMS. 

On s'interesse ici a une forme de calculs dite par blocs, mettant en 
oeuvre la transformation de Fourier ou la transformation de Hartley, cette 
forme necessitant un volume de calculs nettement inferieur ^ la forme 
temporelle dassique. 
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La figure 23 est un schema typique de realisation d*un annuleur 
d'echo sous forme par blocs. Ce dispositif presente des modules C1 et C2 
d'application respectivement d'une contrainte de calcul de correlation et 
d'une contrainte de convolution lineaire. On notera que ces modules CI et 
5 C2 appliquent ces contraintes directement sur un signal temporel et sont 
done chacun encadres par un module de transformee inverse et un module 
de transformee directe. 

La figure 24 montre un schema de realisation d'un annuleur 
d'echo sous forme par blocs selon invention. 
10 Dans ce dispositif. les modules CI et C2 ainsi que les modules de 

transformees inverse et directe les encadrant. ont ete remplaces a chaque 
fois par un module de lissage ne comportant pas de calcul de transformee 
inverse ou directe. 

Le tableau suivant detaille le volume de calculs a effectuer pour 
15 traiter un bloc d'echantillons reels dans le cas d'un module de contrainte 
classique a transformees de taille T et dans le cas d'un filtre de lissage 
frequentiel de 7 coefficients utilisant les reductions de calcul proposees 
precedemment : 





forme classique 


Operation 


forme modifiee suivant 






effectuee 


I'invention 


CI 


2Tlog2T 


(X) 


3T 




2Tlog2T 


(+) 


4T 


C2 


4Tlog2T+T 


(X) 


3T 




4Tlog2T 


(+) 


4T 



2 0 Ainsi, si Ton considere une transformee de taille T = 1024. il vient : 





forme classique 


Operation 


forme modifiee suivant 






effectuee 


I'invention 


C1 


20x1024 


(x) 


3 X 1024 




20x1024 


(+) 


4 X 1024 



C2 


41 X 1024 


(X) 


3 X 1024 




40 X 1024 


(+) 


4 X 1024 



L'emploi du dispositif objet de I'invention represente done un gain 
de 109568 operations arithmetiques. toutes operations confondues (le 
nombre d'op6rations etant ramene a 14336). 

II est bien connu egalement que dans les applications de 
transmission de la parole, le bruit ambiant capte par le microphone du 
terminal peut etre une cause de gene pour les usagers distants (le niveau de 
bruit peut etre tres eleve lorsque la conversation telephonique a lieu depuis 
un poste radiotelephonique mains-libres installe dans un vehicule). 

De meme, pour un systeme de transmission installe dans une 
salle. il est bien connu qu'en absence d'un dispositif annuleur d'echo suivant 
le principe de la figure 22 ou meme en presence d'un dispositif annuleur 
d'echo dont le filtre adaptatif a une longueur sensiblement plus faible que la 
reponse acoustique de la salle, Techo present sur la voie de parole emise 
peut etre une cause de gene importante pour ces memes usagers distants. 

Dans I'etat de I'art classique, des techniques basees sur l'emploi 
d'un attenuateur variable (dites a variation de gain) sont employees pour 
diminuer la gene due a I'echo. Dans I'etat de I'art plus avance. des filtres 
variables ajustes suivant les caract§ristiques des differents signaux (dits 
filtres optimaux) sont employes dans ce but et offrent une qualite subjective 
meilleure, Ces memes filtres variables peuvent etre employes pour reduire 
efficacement le bruit present sur la parole transmise. 

La figure 25 represente un dispositif base a la fois sur Temploi 
d'un annuleur d'§cho et sur Temploi d'un filtre optimal pour effectuer la 
reduction d'une perturbation composee de bruit et/ou d'echo. 

Le module d'annulation d'echo de ce dispositif selon Tinvention 
comporte un haut-parleur 100. un micro 200, un module de soustraction 300 
et un ensemble de filtrage 400. 



Cet ensemble de filtrage 400 regoit a une premiere entree un 
signal x(t) arrivant egalement sur le haut-parleur 100, ce module de filtrage 
recevant Egalement par une seconde entree un signal qui est obtenu en 
sortie du module de soustraction 300, par soustraction sur un signal de micro 
y(t) transform^ dans le domaine frequentiel d'un signal obtenu par application 
du module de filtrage sur le signal de haut-parleur x(t). 

Le module d'annulation d'6cho comporte deux branches se 
rejoignant en un multiplicateur 450. Une premiere de ces branches comporte 
un module 410 de transformee de Fourier directe du signal de haut-parleur 
x(t), et un module de conjugaison 420 de ce signal transforme. La deuxieme 
de ces branches comporte un module de contrainte CI realisant une 
contrainte sur un signal Y' pour calculer un vecteur erreur. un module 440 de 
multiplication par une matrice de normalisation en sortie de ce module de 
contrainte CI. Le module multiplicateur 450 realise la multiplication entre le 
vecteur de haut parleur transforme et conjugue et le vecteur erreur multiplie 
par la matrice de normalisation. 

Ce module 450 de multiplication fournit un vecteur frequentiel a un 
module de contrainte C2. En sortie de ce module de contrainte C2 est place 
un module 470 de mise a jour des coefficients d*un filtre adaptatif de 
I'annuleur d'echo. Le filtre ainsi mis a jour est applique a la transformee du 
signal de haut-parleur x(t) au niveau d'un module multiplicateur 480 dispose 
a la sortie du dispositif de mise a jour 470 et en amont du soustracteur 300. 

Le dispositif selon I'invention presente un avantage fondamental 
sur les dispositifs connus, qui est que le signal frequentiel auquel le filtre est 
applique au niveau du module 480 n'est pas transforme immediatement dans 
le domaine temporel pour etre soustrait au signal temporel d'entree de micro 
y(t), mais est conserve dans le domaine frequentiel. Ainsi, dans ce dispositif, 
c'est le signal de micro y(t) qui est transforme dans le domaine frequentiel 
pour pouvoir soustraire a ce signal de micro transforme dans le domaine 
frequentiel le signal d'echo estime, obtenu en sortie du multiplicateur 480. qui 
est conserve dans le domaine frequentiel. Cette operation de soustraction 



dans le domaine frequenliel est realisee au niveau du soustracteur 300 qui 
est dispose entre le module 480 d'application du fillre adaptatif d'annulation 
d'echo et un module de Iransformee de Fourier directe de sortie de micro 
495. 

Au sein de ce dispositif. ce n'est pas un signal d'erreur tempore! 
correspondant au signal de micro auquel est soustrait le signal d'echo 
estime. qui est transmis au module de contrainte C1. On transmet au module 
de contrainte C1 un vecteur frequentiel d'erreur sans retour dans le domaine 
temporel entre le calcul de ce vecteur d'erreur et le module de contrainte Cl. 
Ce vecteur frequentiel d'erreur est egal a la difference entre la transformee 
frequenlielle du signal de haul parleur auquel est applique le fillre adaptatif 
470, et le signal de micro transforme dans le domaine frequentiel. Ainsi. le 
module de contrainte CI ne presente pas comme dans le cas des dispositifs 
connus, deux sous modules, mais trois sous modules qui sont 
successivement un premier sous-module 432 de transformee de Fourier 
inverse dans le domaine temporel. un second sous-module 435 d'annulation 
d'une partie initiate du signal temporel obtenu en sortie du module 432. et un 
troisieme sous module 438 realisant une transformation de Fourier directe du 
signal temporel ainsi tronque dans le domaine frequentiel. 

Ainsi. le module de contrainte CI realise sur le vecteur d'erreur 
frequentielle qui lui est transmis, un traitement Equivalent a une troncature de 
la transformee temporelle de ce signal d'erreur frequentielle. 

Le signal d'erreur frequentielle Y' est done utilis6 par I'ensemble 
forme des modules 430 de contrainte C1, 460 de contrainte C2. 450 de 
multiplication et 420 de conjugaison, pour calculer. avec un nouveau bloc 
d'echantillons x((t+1 )B). un nouveau jeu de coefficients du filtre adaptatif qui 
sera applique a ce nouveau bloc de donnees de haut parleur x((t+1 )B) pour 
former un nouveau vecteur Y*. 

Le signal Y' obtenu en sortie du soustracteur 300 n'est pas 
seulement foumi au module de contrainte C1. il est egalement transmis a la 
seconde partie du dispositif formant un reducteur de perturbations. 
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Celui-ci est de forme classique en lui-meme, Le signal Y'*obtenu 
en sortie de I'annuleur d'echo est fourni a I'entree d'une boucle 500 decrivant 
successivement deux modules 520 et 530. Le module 520 realise un calcul 
des coefficients du filtre reducteur de perturbations ^ partir du signal 
5 frequentiel Y'. Le vecteur frequentiel de sortie de ce premier module 520 est 
fourni ^ un module 530 realisant une contrainte C3 sur ce vecteur frequentiel. 
Le module 530 se decompose en trois sous modules, un premier sous 
module 532 realisant une transformation de Fourier inverse du vecteur 
frequentiel de sortie du module 520 de fagon a fournir le signal temporel 

10 correspondent a ce vecteur. Ce signal temporel est ensuite soumis a un sous 
module 536 qui annule une partie de ce vecteur et plus precisement une 
partie finale de celui-ci, puis le vecteur temporel tronque subit une 
transformee de Fourier directe au niveau d'un troisieme sous module 538. La 
boucle 500 fournit ainsi un filtre reducteur de perturbations dont les 

15 coefficients frequentiels sont tels que leur transformee de Fourier inverse a 
un nombre de coefficients non nuls reduit. 

Le vecteur frequentiel obtenu en sortie du second module 530 est 
applique au signal frequentiel d'erreur Y' au niveau d'un module de 
multiplication 600. On obtient ainsi un signal frequentiel Z' correspondant a la 

2 0 transformee de Fourier d*un signal temporel dont I'echo a ete annule et dont 
des perturbations, ici le bruit, ont ete reduites. 

Le vecteur frequentiel Z' ainsi obtenu est ensuite transforme dans 
le domaine temporel par un module de sortie 700 realisant une transformee 
de Fourier inverse sur ce signal Z'. On obtient alors un bloc d'echantillons 

2 5 temporel dont I'echo a 6te annule et dont les perturbations ont ete reduites. 

On notera egalement que Ton dispose en sortie de ce dispositif un 
module 800 realisant un passage d'une cadence bloc a une cadence 
echantillon. Le message sonore obtenu est ensuite transmis a un 
correspondant qui regoit done le signal sonore present devant le micro 200 a 

3 0 {'instant considere, duquel ont ete supprimes I'echo et le bruit. 
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Ainsi ce dispositif permet de supprimer un module de Iransformee 
de Fourier directs qui existait classiquement entre un module de calcul 
d'erreur et un module de reduction de perturbations lorsque deux tels 
modules etaient coupl6s. Le signal fr6quentiel est done fourni directement 
^ la boucle 500 et notamment au premier module de cette boucle realisant 
un calcul du filtre reducteur de perturbations. 

On obtient ainsi en sortie du module de contrainte 530 un jeu de 
coefficients du filtre reducteur de perturbations adapte en fonction du signal 
Y\ 

Le dispositif qui vient d'etre deait presents un reducteur de bruit. 
Toutefois ce reducteur de bruit peut etre remplace par un module de 
reduction d'echo. un reducteur diminuant un melange d'echo et de bruit, ou 
plus generalement tout reducteur de perturbations. 

On notera la presence de trois modules de contraintes CI. C2 et 
C3 consistant chacun en un module d'annulation d'une partie d'un signal 
temporel encadre par deux transformees, une inverse et une directe. La 
contrainte C3 peut etre vue comme une contrainle de convolution lineaire. 
Selon rinvention. on peut considerer C3 comme un moyen pour lisser la 
reponse en frequence du filtre et controler la selectivite de ce filtre en limitant 
son support temporel et en annulant done autant d'echantillons que 
necessaire pour effectuer ensuite la transformee de taille T. 

Conformement a I'invention, on remplace chacun des trois 
modules de contrainte CI, C2. C3 par un module de lissage frequentiel. 
Selon le but recherche, la contrainte C3 peut etre remplacee par un module 
de lissage frequentiel ou un module d'interpolation faisant objet de 
I'invention. Le volume des calculs effectues dans le dispositif de rinvention 
est nettement plus faible que dans le cas du dispositif classique. 

Ainsi, sur la figure 25. la realisation d'un annuleur d'echo combine 
a une reduction de perturbations par filtrage optimal necessite I'emploi de 9 
transformees. Le dispositif de la figure 26 ne differe de celui qui vient d'etre 
decrit en reference a la figure 25 qu'en ce que les trois modules de contrainte 
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CI, C2 et C3 ont ete remplaces par des modules de lissage utilisant une 
convolution avec une fonction U telle que celles qui ont ete decrites 
precedemment. 

Comme il est montre sur la figure 26, un syst^me correspondant 
selon rinventlon necessite seulement 3 transformees et 3 modules de 
lissage, les trois modules de lissage remplagant chacun un ensemble de 
deux transformees et d'une contrainte temporelle. 

Le tableau suivant detaille le volume de calculs a effectuer pour 
realiser les contraintes CI, C2 et C3 dans le cas d'une transformee de taille 
T (figure 25) et dans le cas d'un filtre de lissage frequentiel de 7 
coefficients (figure 26) : 





Forme classique 


Operation 


Forme modifiee suivant 






effectuee 


("invention 


CI 


4Tlog2T + T 


(X) 


4T 




4Tlog2T 


(+) 


. 7T 




4Tlog2T + T 


(X) 


4T 




4Tlog2T 


(+) 


77 


C3 


4Tlog2T + T 


(X) 


AT 




4Tlog2T 


(+) 


7T 



Ainsi, si Ton considere une transformee de taille T = 1024, il vient : 





Forme classique 


Operation i Forme modifiee suivant 






effectuee 


j I'invention 


CI 


41 X 1024 


(X) 


4 X 1024 




40 X 1 024 


(+) 


7 X 1024 


C2 


41 X 1024 


(X) 


4 X 1024 




40 X 1 024 


(+) 


7 X 1024 


C3 


41 X 1024 


(X) 


4 X 1024 




40 X 1024 


(+) 


7 X 1024 



L'emploi d'un dispositif conforme a I'invention represente done un 
gain de 215040 operations anthmetiques. toutes operations confondues. (Le 
nombre d'operations etant ramene a 33792). 

II permet de faire directement dans le domaine transforme. c'est a 
dire le domaine frequentiel. Tequivalent de la mise a zero dun nombre donne 
de coefficients temporels (lissage) ou Tequivalent de I'ajout d'un nombre 
donne de zeros dans le domaine temporel (interpolation), condition 
necessaire pour ne pas generer de convolution/correlation circulaire lorsque 
Ton fait du filtrage/correlalion dans le domaine transforme, comme cela est 
expose en detail dans le document [6]. 

Ainsi. un dispositif conforme a celui de la figure 25 est 
particulierement adapte a une reduction importante de calculs en utilisant 
une fonction U telle que celles decrites precedemment. Un tel dispositif. muni 
de tels modules de lissage, permet d'effectuer I'annulation d'echo et la 
reduction de perturbations totalement dans le domaine frequentiel, sans 
retour dans le domaine temporel couteux en calculs. 

Les differents dispositifs decrits peuvent etre realises au moyen de 
la technique de Tetat de I'art, notamment au moyen de processeurs de signal 
en temps reel (DSP) typiquement a arithmetique flottante, tel que ceux 
couramment utilises dans les applications audio et de telecommunications 
(par exemple TMS320C3X ou TMS320C4X de TEXAS INSTRUMENTS) ; 
des processeurs a arithmetique fixe peuvent aussi etre utilises, en respectant 
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les precautions usuelles de cadrage approprie des donnees dans les 
traitements. II est a noter que Tinvention ne presents pas uniquement des 
avantages en complexity, mais elle rend aussi rimpl6mentation plus facile 
(pas de permutation « scrambling », moins de probl6mes de cadrage pour 
calcul en virgule fixe, programme plus simple, adressage plus simple,...). 

Un dispositif selon Tinvention est done du plus grand interfet des 
que Ton fait du filtrage variable dans le temps et que Ton travaille dans le 
domaine transforme. 

Ceci englobe done non seulement le filtrage adaptatif utilise en 
annulation d'echo mais notamment, de fa?on plus generale, tout filtrage 
adaptatif dont les utilisations sont multiples : debruitage, separation de 
signaux. etc. 

Un contexte typique d'application d'un dispositif de I'invention est 
decrit dans [7], ou un dispositif classique a deux transformees encadrant un 
module de fenetrage est utilise pour attribuer au filtre des caracteristiques 
selectives en modifiant la taille de la reponse impulsionnelle de celui-ci, 

Dans certains modes de realisation decrits ici, on convolue avec U 
une fonction X qui est la Transformee de Fourier discrete d*un signal 
temporel discret. Bien entendu. une fonction U telle que celles decrites ci- 
dessus s'applique de la meme fagon a des signaux frequentiels 
correspondant a une transformee de Fourier discretisee d*un signal temporel 
continu. 
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REVENDICATIONS 

1. Dispositif de traitement de donnees numeriques comportant un 
module (M2, M3, M'a) realisant sur un vecteur de donn6es du domaine 
frequentiel Z(k). ou k varie de 0 ^ N-1, une convolution avec una fonction U, 
convolution qui correspond ^ une annulation dans la domaine tempore! 
d'echantillons de la transformee inverse de Z(k), caracterise en ce que la 



fonction U est de la forme : Uik) = sine 



— ^je ' \P{k) . ou koest un 



10 entier constant et P(k) est une fenetre de ponderation symetrique autour de 
ko. 

2. Dispositif selon la revendication precedente, caracterise en ce 
que ko est egal a zero. 

3. Dispositif selon I'une des revendications precedentes, caracterise 
15 en ce qu'il comprend un module (M3') recevant un vecteur frequentiel (H), 

aple a inserer entre deux coefficients de ce vecteur (H) a chaque fois un 
coefficient supplementaire de fa^on a fournir un vecteur frequentiel (H') de 
longueur augmentee. 

4. Dispositif selon la revendication precedente, caracterise en ce 
20 que les coefficients supplementaires sont des zeros. 

5. Dispositif selon i'une des revendications 3 ou 4, caracterise en ce 
que, le vecteur frequentiel de longueur augment6e (H*) etant indice de 0 a 
2N-1, les coefficients inseres sont les coefficients d'indices impairs. 

6. Dispositif selon I'une des revendications 3 a 5. caracterisd en ce 
2 5 que le module (M2, M3, M"3) realisant une convolution avec U est place en 

aval du module d'insertion (M'3) et en ce que Z est le vecteur frequentiel de 
longueur augmentee (H'). 

7. Dispositif selon I'une quelconque des revendications precedentes. 
caracterise en ce qu'il comporte un filtre (H) en amont du module 



(M2.M3.M"3) de convolution avec U, et en ce que Z est la fonction de 
Iransfert de ce filtre (H). 

8. Dispositif selon la revendicalion precedente. caracterise en ce qu'il 
comporte des moyens de calcul (H) des coefficients du filtre (H) a partir d'un 
signal d'entree (X.S1) du dispositif. 

9. Dispositif selon Tune quelconque des revendications prec6dentes. 
caracterise en ce qu'il comporte un module (Ml) apte a realiser un 
transformee dans le domaine frequentiel d'un signal temporel d'entree (SI), 
ce module (Ml) etant en amont du module {M2.M3.M"3) de convolution avec 
U, et en ce que Z est cette transformee frequentielle (X). eventuellement 
completee par Tintermediaire d'un module d'insertion (M's) conforme a Tune 
des revendication 3 a 5. 

10. Dispositif selon les revendications 6, 7 et 9 en combinaison. 
caracterise en ce qu'il comporte un premier module (M2) appliquant une 
convolution avec une premiere fonction de la forme 



t;(/c) = sinc|^— j-^je ^ ^ rP{k) sur une transformee frequentielle (X) 

d'un signal d'entree (S1) eventuellement complete, un second module de 
convolution (M3) appliquant une convolution avec une deuxieme fonction de 



la forme L/(/c) = sin c ° [e ^ ^ \P{k) sur la reponse frequentielle (H) 



d'un filtre (H) eventuellement completee. les vecteurs de sortie (X'.S3) de ces 
deux modules (M2.M3,M"3) ayant un meme nombre de coefficients, et en ce 
que le dispositif comporte en sortie un module (M4) apte a multiplier les 
coefficients de ces deux vecteurs de sortie (X',S3). 

11. Dispositif selon Tune des revendications precedentes en 
combinaison avec Tune des revendications 3 a 5, caracterise en ce que le 
module (M2,M3,M"3) de convolution avec U fournit un vecteur de sortie 
(H\S3) de meme longueur que le vecteur complete Z, en conservant dans le 
vecteur de sortie (H\ S3) les coefficients du vecteur Z qui etaient presents 
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avant insertion, et les autres coefficients du vecteur de sortie (H\S3) etant 
obtenus par convolution de Z et U. 

12. Disposilif selon Tune des revendications pr§c§dentes. 
caracterise en ce que le module (M2,M3,M' 3) de convolution avec U fournit 
5 en sortie un vecteur B(0,...,N-1) qui est tel que quel que soit k. un coefficient 
d'indice k de B est egal a un produit de convolution entre Z et U qui est tel 
que le coefficient d'indice k de Z est multiplie dans ce produit de convolution 
avec le coefficient d'indice ko de U pour lequel la fonction sine a pour 
argument 0. 

10 13. Dispositif selon Tune quelconque des revendications 

precedentes, caracterise en ce que le filtre U prend des valeurs non nuiles 
sur un intervalle de valeurs de k qui est symetrique par rapport a la valeur ko 
pour laquelle U atteint son maximum de module. 

14. Dispositif selon Tune des revendications precedentes. caracterise 
15 en ce que la fonction U comporte un nombre de coefficients Lu impair, et en 

ce que U peut s'ecrire : U{k) = sin c + " k ^\P(k) 

\ A 2J 

15. Dispositif selon Tune quelconque des revendications 
precedentes en combinaison avec Tune des revendications 9 ou 10. 
caracterise en ce que la transformee est une transformee de Fourier 

20 discrete. 

16. Dispositif selon Tune des revendications precedentes, 
caracterise en ce que la fenetre de ponderation est une fenetre de Kaiser de 
coefficient 1 ,5. 

17. Dispositif selon I'une des revendications precedentes, caracterise 
2 5 en ce qu'il constitue un annuleur d'echo. 

18. Dispositif selon Tune des revendications precedentes, caracterise 
en ce qu'il constitue un reducteur de bruit. 

19. Dispositif selon Tune des revendications precedentes, caracterise 
en ce que a = 1 . 
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20. Dispositif selon Tune des revendicalions precedentes, caracterise 
en ce que a = -1 . 

21. Dispositif selon Tune des revendications prec^dentes, caracterise 
en ce qu'il comprend un haul-parleur (100). un micro (200). un annuleur 
d^echo (420,430.440.450) et un reducteur de perturbations (500). I'annuleur 
d'echo comportant un filtre adaptif (470) et un module de soustraction (300) 
fournissant I'erreur (V) entre un signal provenant du micro (200) et un signal 
obtenu par application du filtre adaptatif (460) a un signal de haut-parleur 
(100), le filtre adaptatif (460) adaptant ses coefficients en fonction de cette 
erreur, et le dispositif comportant des moyens (495) aptes a transformer le 
signal de micro dans le domaine frequentiel en amont du module de 
soustraction (300) de fafon a ce que la soustraction soit effectuee dans le 
domaine frequentiel. 

22. Dispositif selon la revendication precedente, caracterise en ce 
qu'il presente des moyens (430.440) pour transmettre le resultat de cette 
soustraction frequentielle au filtre adaptatif (470) de I'annuleur d'echo. 

23. Dispositif selon la revendication 21. caracterise en ce que le 
reducteur de perturbations (500) est place en aval du module soustracteur 
(300) et est applique dans le domaine frequentiel au resultat de la 
soustraction. 

24. Dispositif selon la revendication precedente, caracterise en ce 
que le reducteur de perturbations (500) comporte un filtre adaptatif (520), 
apte a recalculer ses coefficients en fonction d'un signal d'entree frequentiel 
(V) du reducteur de perturbations (500). 

25. Dispositif selon la revendication precedente, caracterise en ce 
que le reducteur de perturbations (500) est place pour recevoir en tant que 
(edit signal d'entree frequentiel du reducteur de perturbations le signal (Y") 
frequentiel de sortie du module soustracteur (300). 

26. Dispositif selon Tune des revendications 24 ou 25. caracterise en 
ce que le reducteur de perturbations (500) forme une boucle recevant en 
entree le signal frequentiel de sortie (Y') du soustracteur (300). et appliquant 
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10 



en sortie une multiplication des coefficients adaptes de son filtre adaptatif sur 
le signal frequentiel de sortie (Y') du soustracteur (300), 

27. Dispositif selon Tune des revendications 24 ^ 26. caracterise en 
ce que le m6me signal frdquentiel (Y') est utilise comma erreur pour adapter 
le filtre adaptatif (470) de I'annuleur d'6cho et est multipli6 par les coefficients 
du filtre adaptatif (520) du r§ducteur de perturbations (500). 

28. Dispositif selon Tune des revendications 21 S 27, caracterise en 
ce qu'aucun module de transformee n'est place entre le module soustracteur 
(300) et le reducteur de perturbations (500). 
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